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ABSTRACT A simple analytic calculation is presented for the secu- 
lar evolution of the eccentricities in a hierarchical triple system such as the 
one containing the millisecond pulsar PSR 1620 — 26 in M4. If the second 
companion of PSR 1620 — 26 is of stellar mass (m2 ^ 0.1 Mq), an eccen- 
tricity as large as that observed today for the inner binary (ei ~ 0.03) 
could very well have been induced by secular perturbations. In contrast, 
such a large eccentricity cannot be induced by a second companion of 
planetary mass. 



INTRODUCTION 

The millisecond pulsar PSR B1620 — 26 in the globular cluster M4 has a low- 
mass companion (mass mi ~ 0.3 Mq for a pulsar mass rrip = 1.35 M©) in a 
nearly circular orbit of period Pi = 0.524 yr (Lyne et al. 1988; McKenna & 
Lyne 1988). The eccentricity ei = 0.0253, although small, is several orders 
of magnitude larger than observed in most other low-mass binary millisecond 
pulsars (Thorsett, Arzoumanian, & Taylor 1993, hereafter TAT). In addition, 
the timing data reveal a very large value of the pulse period second derivative 
P = -2.3 X 10-"s-i (Backer 1993), indicating the presence of an additional 
source of acceleration. It is very unlikely that this could be caused by a nearby 
passing star or by the mean gravitational field of M4. Instead, the most likely 
explanation is that the pulsar has a second, more distant orbital companion 
(Backer 1993; Backer, Foster, & Sallmen 1993; TAT). Such hierarchical triple 
systems are expected to be produced quite easily in globular clusters through 
dynamical interactions between binaries (Mikkola 1984; Hut 1992; Rasio, Hut, 
& McMillan 1994). 

The present timing data are consistent with a second companion mass any- 
where in the range 10"'^ Mq ^ m2 ^ 1 Mq with a corresponding orbital pe- 
riod 10 yr ^ P2 ^ 10'^ yr (see Michel 1994, and these proceedings). Thus, in 
particular, there is a possibility that the second companion might be a Jupiter- 
like planet (Backer 1993; Sigurdsson 1993; TAT), although it could also be an- 
other star (main sequence, white dwarf, or even another neutron star). Planets 
have been detected in orbit around at least one other millisecond pulsar, PSR 
B1257+12 (Wolszczan & Frail 1992; Wolszczan 1994). 

If PSR B1620 — 26 is indeed in a triple system, then small perturbations 
of the (inner) binary pulsar's orbit should be induced by the presence of the 
second (more distant) companion. Rasio (1994) has argued that the unusually 



large eccentricity of the binary pulsar could be explained naturally as arising 
from these perturbations, but only if the second companion is a star and not a 
planet. A completely analytic derivation of the induced eccentricity is given here 
for various limiting regimes of interest. Other treatments of secular evolution in 
triple systems, using very different approaches, are given by Mazeh & Shaham 
(1979) and Bailyn (1987). 



SECULAR PERTURBATION THEORY 



Assume for simplicity that the orbital eccentricities and the relative inclina- 
tion are low. Standard secular perturbation theory of celestial mechanics (e.g., 
Brouwer & Clemence 1961; see Dermott & Nicholson 1986 for a useful summary) 
can then be used to calculate the subsequent evolution of the system over a time 
t ^ P2 > Pi- The general solution for the eccentricities ej(^t) and longitudes of 
pericenter u)j(^t) can be written 

: e'^+hm{g+t + (3+) + e'fhm{g_t + (3_), (1) 

: e'^+\os{g+t + (3+) + e'f\os{g_t + (3_), (2) 

with j = 1 for the inner orbit and j = 2 for the outer orbit. This solution can 
be interpreted geometrically as the sum of two vectors describing circles in the 
(ej cosujj, ej sin Wj) plane (see Fig. IV of Malhotra, these proceedings). Here g± 

and e^-^^ are the eigenvalues and corresponding eigenvector components of the 
2x2 secular perturbation matrix with components 

All = ^TCi/X2a26^y2(") = o-?", (3) 

A12 = -^TCi/xaa^fe^^^^^C") = -(^QctP, (4) 

A21 = -\n2lJ,iab''^j'^{a) = -avfi, (5) 

A22 = \n2^llab^^j^{a) = au, (6) 

where a = ai/a2 < 1 is the ratio of semimajor axes, Uj = 27r/Pj, /Xj = mj/mp, 
and a = [l/4:)niiJ,iab^^2{^)- The dimensionless functions b^j^\a) are Laplace 
coefficients (Brouwer & Clemence 1961). Following the notations of Rasio et 
al. (1992), we have introduced the quantities q = /X2//X1, 1/ = n2/ni, and /3 = 

^3/2/^3/2' '^^^ eigenvalues of the matrix Aij can then be written explicitly 

5± = I + ± [(?" - i^f + ^0iuf3^\ '^'1 , (7) 

and the eigenvector components are given by 

.(+) 



up to a normalization. This normalization, as well as the phases /3+ and /3_ 
appearing in equations (1) and (2), must be determined from an initial condition 
(i.e., values of ej and u)j at t = 0). 

Here we assume that the (inner) binary pulsar had a very low eccentricity 
ei(0) ~ initially, i.e., at the time it acquired its second companion. One can 
then set arbitrarily wi(0) = ^2(0) = 0. It is straightforward to show that the 
corresponding solution of equations (1) and (2) for t = implies /3+ = /3_ = 
and eigenvectors with 

(+) _ e(-)e2(0) (-) _ -62(0) . . 

Solving equations (1) and (2) in general for ei(f) then gives 

ei(f) = J^{q, a) 62(0) [1 - cos{gt) ]^/\ (11) 
where g = — and 

e(-) 



1 - e(+)e(- 



(12) 



(cf. Rasio 1994). Thus the eccentricity of the inner binary returns periodically 
to its initial value ei(0) ~ 0, with a maximum amplitude ei^^ax = \/2-^e2(0) 
proportional to the eccentricity of the outer orbit. 

The analytic solution given by equations (7)-(12) is valid only for low ec- 
centricities and low relative inclination. In general, numerical integrations of the 
three-body problem must be used (Rasio 1994). One finds that the form of the 
solution remains generally similar to that of equation (11), with ei(f) returning 
to zero periodically. However, both the period and the amplitude of the solution 
can deviate significantly from their analytic values when 62 or the inclination is 
very large. 



LIMITING SOLUTIONS FOR a < 1 

In a strongly hierarchical triple system one can assume that a <C 1. The relevant 
Laplace coefficients can then be expanded as b^^^ = 3a + 0{a'^) and b^^^ = 
(15/4)a2 + C)(a^) (Brouwer & Clemence 1961), so that /3 ~ (5/4)a. Using 
jy ~ a'^/^ we get from equation (7) 

5± ~ I i^qa + ± [q'a' + - 2qa'^'] '^'} . (13) 

To proceed further, we must distinguish between two different limits, corre- 
sponding to q ^ a^l'^ and q <C a-^/^. We refer to these limits as the "stellar" 
and "planetary" cases, respectively. 

In the stellar case one finds that the eigenvalues given by equation (13) can 
be written 

5+ ~ aqa, g_ ~ aa^l'^ < 5+, (14) 



and the corresponding eigenvectors, calculated using equations (8) and (9), have 
e(+)^_Aa3/2 e(-)~-a. (15) 

Equation (12) then gives 

5 3 
^{q,Oi)c^ ^-^a, 5 ~ 5+ ~ -TCi/X2a^. (16) 

Remarkably, the maximum induced eccentricity is independent of the second 
companion mass in this limit. However, this maximum eccentricity is reached 
after a time inversely proportional to /X2. 

For the planetary case the eigenvalues are 

5+ ~ cra^/^ 5f_ ~ CT^a < 5+, (17) 

and the eigenvectors have 

~ _ A a-i/2 g(-) ^ 1 a-i (is) 

giving 

5 3 

^{q,Oi)^ ^^qa'^^^, g 2^ g+ 2^ -uifj.ia'^/^. (19) 

Here the period of the eccentricity variations is independent of /X2, but the max- 
imum amplitude decreases like the mass ratio q. 



APPLICATION TO THE PSR 1620-26 TRIPLE 

For PSR 1620 — 26 we have a ^ 10~^ and q ~ 10"'^ for a Jupiter-size second 
companion (planetary case) or ~ 1 for a stellar-mass second companion. Using 
equations (14)-(19) we obtain for the maximum induced eccentricity ei^^ax = 
V2Te2, 

( 3 X 10-^62 (a/10-2) stellar case , , 

ei,max - I 3 X 10-4 62 {q / 10-^){a / 10-^^/^ planetary case ^ 

Clearly, if the present eccentricity ei = 0.025 has been induced by secular per- 
turbations, the second companion must be of stellar mass (Rasio 1994). 

The period of the eccentricity variations is Pg^ = 2Tr/g, and we find in the 
two limiting cases 

^ f 7 X lO^yr/x^^ (a/10-2)-^ stellar case , . 

~ [ 3 X 10^yr(a/10-2)-^/2 planetary case ^ 

This is certainly much shorter than the age of the pulsar, tp ^ 10® yr (TAT). 
It is also shorter than the mean time between close encounters with passing 
stars in the core of M4 (density p = lO'^Pi Mq pc"'^ and velocity dispersion 
a = 5c75kms--^), which is tc ~ lO^yip^ C75(a2/10 au)"-^. 
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